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A mixture of light and heavy spin-polarized fermionic atoms in an optical lattice is 
considered. Tunneling of the heavy atoms is neglected such that they are only subject to 
thermal fluctuations. This results in a complex interplay between light and heavy atoms 
caused by quantum tunneling of the light atoms. The distribution of the heavy atoms is 
studied. It can be described by an Ising-like distribution with a first-order transition from 
homogeneous to staggered order. The latter is caused by an effective nonlocal interaction 
due to quantum tunneling of the light atoms. A second-order transition is also possible 
between an ordered and a disordered phase of the heavy atoms. 

1. Introduction: Phase separation and the physics of correlated disorder 

Mixtures of different species of cold atoms open an interesting field of many-body 
physics, where the competition between the atoms may lead to new quantum states. This 
is the case for low-temperature properties (i.e. for the ground state) as well as in the 
presence of thermal fluctuations. For the latter case we expect phase separation, where 
different phases can coexist. This phenomenon has been intensively discussed in the solid- 
state community in the context of complex materials pQ, and it is believed to be a result 
of the competition of different (e.g. spin and orbital) degrees of freedom j2]. A similar 
situation may appear in atomic mixtures due to correlations caused by the competition 
of the different atomic species. Two or more phases can coexist if there is a first-order 
phase transition between the individual phases. In other words, thermal fluctuations 
lead to metastable states that decay only on long-time scales. Such a phenomenon is 
characterized by a broken translational invariance. In terms of an atomic mixture in an 
optical lattice the atoms are arranged in clusters, consisting of one of the phases, with a 
characteristic length scale relative to the optical lattice constant. This implies a specific 
type of correlated disorder that has interesting dynamical properties of the mixtures, 
including the localization of atoms. 

The physics of disorder in statistical or solid-state physics is associated with random 
scattering of light particles by a disordered (i.e. a randomly produced) array of slow 
particles. This disordered array is static and does not change during the course of the 
considered scattering process. For instance, the scattering of electrons in a crystal with 
randomly distributed impurity atoms is a typical realization of disordered physics or the 
scattering of photons in atmospheric dust. 

A comparable situation is found in an atomic system subject to an optical lattice. The 
periodic structure of the optical lattice, given by counterpropagating lasers, is normally 
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quite robust such that disorder in the sense of impurities cannot easily be achieved. 
However, in contrast to a crystal in solid-state physics, two or more different types of 
atoms can be mixed and brought in an optical lattice, where one type of the atoms is 
relatively light (e.g. 6 Li) and the other type(s) is (are) heavy (e.g. 40 K, 23 Na [3], 87 Rb). 
The different masses lead to different dynamical properties. In particular, heavy atoms 
behave almost like static degrees of freedom, in comparison to the light atoms. Therefore, 
heavy atoms play the role of the impurity atoms in the crystal and the light atoms the 
role of the electrons. Thus disorder physics can be studied in mixtures of atoms with 
different masses which are available in experiments The main difference is here that 
we shall study the interaction of two species of atoms within the same model. In other 
words, the distribution of disorder (i.e. the distribution of the heavy atoms) is a direct 
consequence of the interaction between light and heavy atoms. 

Disordered systems have attracted a lot of attention because they provide a new class 
of physics, including new phases and new types of phase transitions. The main reason is 
that scattering of quantum particles by a periodic structure is qualitatively different from 
scattering by a disordered structure: The former is governed by Bloch's Theorem [5 J which 
states that the scattered quantum states are not strongly affected by the scattering process 
but keep their wave-like properties of the unscattered states. This changes dramatically 
as soon as the periodicity of the scatterers is disturbed. At first the quantum states loose 
their wave-like character because the physics is controlled by diffusion. And for stronger 
disorder the quantum states can even be localized P • There is a phase transition between 
the diffusion-controlled regime and the localized regime jjj. 

A fundamental concept for a theoretical description of disordered system is that physical 
quantities are averaged with respect to a statistical distribution of the random structures. 
It will be discussed in this article that averaging over thermal fluctuations of heavy atoms 
provides such a type of disorder (or quenched) average of physical quantities, like the 
density or the dynamical Green's function of a light atom. Emphasize is on fermionic 
system, since the appearence of Bose-Einstein condensation should not be included in 
the discussion. Moreover, there is a formal mapping of the heavy fermions to Ising spins 
(correlated binary alloy). This provides the opportunity to discuss the statistics in terms 
of para-, ferro- and antiferromagnetic Ising states. 

The article is organized as follows: In Sect. 2 the model of a mixture of two types of 
atoms is discussed for light fermionic and heavy fermionic or bosonic atoms. The effective 
distribution of the heavy atoms is studied in Sect. 3, where it is shown that for fermionic 
atoms this distribution is related to a classical Ising-spin model. The distribution of 
the Ising spins is investigated in Sect. 4, using a classical limit (Sect. 4.1) and in a 
strong-coupling expansion (Sect. 4.2). The results are discussed in Sect. 5. 

2. Model 

A mixture of two types of atoms is considered. The atoms are subject to thermal 
fluctuations that are treated within a grand-canonical ensemble. Moreover, it is assumed 
that the atoms are in a magnetic trap such that they are spin polarized. There is a local 
repulsive interaction between the two species but no other interaction of the atoms within 
the same species. 
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(c) are creation (annihilation) operators of the light fermionic atoms, /' (/) are 
the corresponding operators of the heavy atoms. The latter can either be fermionic or 
bosonic. This gives the formal mapping 



3 Li — ► 4, c T 40 K ( 23 Na, 87 Rb) — > f r , f r . 



The physics of the mixture of atoms is defined by the (asymmetric Hubbard) Hamiltonian: 



H = -t c E 4<V - tf E flfr> + E [ - fdCr - Vffrfr + Uflf r c\i 
(r,r') {r,T') r 



An effective interaction among fermionic species is controlled by the (repulsive) Pauli 
principle. In the case of heavy bosonic atoms the interaction is neglected. 

The tunneling rate decreases exponentially with the square root of the mass of the 
particle. If the / atoms are heavy, the related tunneling rate is approximated by if ~ 0. 
This limiting model is known as the Falicov-Kimball model |8|9|lUj . 

A grand-canonical ensemble of atoms in an optical lattice is a situation in which the 
optical lattice can exchange atoms with a surrounding atomic cloud. This is given at the 
temperature 1/(3 by the partition function 

Z = Tre-? 11 . 

The real-time Green's function of light atoms describes the motion of a single light atom, 
interacting with the atomic mixture. For the motion from lattice site r' to r during the 
time t it is defined as 



G c (r,it;r',0) = — Tr 



'-PH itHJ -itH 



(2) 



The density of heavy atoms n f jT and of light atoms be obtained from the Green's 

function in the special case t — 0, r' = r: 



n 



:Tr 



-m ft 



frfr 



;Tr 



(3) 



3. Ising representation of heavy fermionic atoms 

The quantities of physical interest (density, Green's function etc.) are of the form of a 
trace with respect to c and / atoms: 

-0H itHJ -OH 



(4) 



If if = this can also be expressed as 
E Tr e (e 



0H{{nr}) c jtH({nr})J p -itH({n r }) 



c r ,e 



{n r } 



since all operators under the trace are diagonal with respect to the number states |{n r }) 
of the / atoms. Therefore, the Hamiltonian depends only on the integer numbers {n r } 
and the c operators as 

H({n r }) = -i E c l°r' + E [nA-^f + Uc\.c r ) - 

(r,r>) r 



= -fif^2n r + H c ({n r }), (5) 
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where t c = t has been used. This means that the / density fluctuations n r have been 
replaced by classical variables: 



frfr 



n, 



which is n r — 0, 1 in the case of fermions and n r = 0, 1, 2, ... in the case of bosons. Thus 
the expression in Eq. (£Q) can be written as a sum over all realizations of n T as 

e^f^r nr Tr c ^ e "^ c({nr}) c r e^ c({nr}) c t ,e" i * Hc({nr}) ) . 

{n r } 

H c ({n r }) is a quadratic form of the c operators. Therefore, it describes independent 
spinless fermions which are scattered by heavy atoms, represented by n r . In particular, 
the density n c , defined in Eq. (JSJ), reads 



{n r } 

Since H c ({n r }) is the Hamiltonian of a noninteracting Fermi gas (it interacts only with 
the density of heavy atoms), the trace Tr c in the partition function can be evaluated as 
a fermion determinant [TT] : 

Z = E e^'ErTr,. (e"^^") = E e ^ ^™ r det[l + e"^], (6) 

{n r } {n r } 

where /i c is an iV x JV matrix for N lattice sites with matrix elements 

where e.,- is a lattice unit vector in direction j. A similar calculation gives for the density 
n c,r = ^ E e^^'detfl + e -^][ e -^(l + e"^)" 1 ]^. (7) 



In Eq. ((7j) there is a non-negative factor 

P({n r }) = -L^^"-det[l + e"^] (8) 

Zj 

which gives J2{n r } P({ n r}) = 1 because of Eq. ©. Thus P({n r }) is a probability distribu- 
tion and the expressions in Eq. (J7|) is a quenched average with respect to this distribution: 



n r 



([e-^l + e-^r'U/- (9) 



The distribution is a realization of correlated disorder and in the case of fermions it is 
a correlated binary alloy. Moreover, the density of heavy atoms is given through the 
distribution P({n r }) as 

n f = (n r ) f = E n r P({n r }). (10) 

{n r } 
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4. Approximations for the distribution of heavy fermionic atoms 

Only fermionic mixtures are considered because their distribution is simpler due to the 
fact that n r has only two values. The latter implies that n r can be expressed by Ising 
spins S r = ±1 as 

S r = 2n r - 1. (11) 

It is possible to discuss the corresponding distribution in terms of magnetic Ising states. 
This will be used in a strong-coupling expansion of the distribution of the heavy atoms 
P({S r }) in powers of t/U. This is an extension of the strong-coupling approximation for 
a fermionic mixture in the symmetric case \i c = fif 



4.1. Classical limit: t — 

The starting point of the strong-coupling expansion (i.e. the unperturbed case) is the 
limit i = 0. This limit is relevant for a mixture in which both atomic species are so heavy 
that tunneling can be neglected completely. The distribution P({n r }) factorizes on the 
lattice in this case and gives 

P({n r }) oc J] e^f nr [l + e ^o-Un r )-\ _ ( 12 ) 

r 

Then the densities of light and heavy atoms are obtained from Q and (|l(Jj) as 

e /3£tc _|_ e /3(M/+Mc-C/) 



X _l_ e /3/ic _|_ e /?Mf _|_ e ^(/"/+Mc-f) 

e /3M/ _|_ e P(Hf+Hc-U) 

X _|_ g/^Mc _|_ e PVf _|_ g/3(M/+Mc-t/) ' 



The densities of the two species are controlled by their chemical potentials. At low 
temperatures (i.e. for (3 ~ oo) they change with fif and fi c in a step-like manner (cf. 
Table 1): The density with the higher chemical potential is always 1 whereas the density 
with the lower chemical potential jumps from to 1 at the value U. This means that the 
repulsive interaction between the different species prevents the fermionic atoms with the 
lower chemical potential to enter the optical lattice unless its chemical potential is strong 
enough to overcome the repulsion. 

For the symmetric case // c = /// = \i the densities are equal: n c = rif. The behavior 
is different in comparison with the asymmetric case fi c ^ fif because of an intermediate 
regime for < \i < U, where both atomic species have to share the optical lattice. As a 
consequence, each density then is 1/2: 



n f = n c = { 1/2 < fi < U . (13) 
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n c 


n f 


> He- 




flc 


< u 





1 




flc 


> u 


1 


1 


ft f < f-c- 




fif 


< u 


1 







ftf 


> u 


1 


1 



Table 1 

Mixture of heavy fermionic atoms in the absence of tunneling: Density of light atoms n c 
and of heavy atoms nf are step-like functions of the corresponding chemical potentials. 
The step occurs when the chemical potential exceeds the repulsive interaction energy U. 



4.2. Strong-coupling (tunneling) expansion: t/U <C 1 

The strong-coupling expansion of the distribution function P({S r }) in powers of t/U 
gives an expansion in terms of the Ising spins |12|13j (cf. Appendix A): 



P{{S r }) oc exp ( Ei S r + S r S r 



m 



1 + e 



f3(tM c -U(l + S r )/2) 



with coefficients 



Ei-- 
and 
E2 - 



P 2 P 



a f3v+l3g 



'I _|_ e Pv-Pg}2 ft _|_ e pu+p g y 



P 2 P 



l3v-/3g 



(3v+(3g 



+ 



smh(/3g) 



Agl + e 2 ^ + 2e? v cosh((3g) 



1 _|_ e Pv-Pg\2 ft _j_ e t3v+0gj2 

These coefficients simplify substantially in the low-temperature regime (j3 ~ 00): 
Ei ~ 0, E2 r 

which gives three different regimes for the distribution: 

He < 



-/3P/4U forO<fi c <U 
otherwise 



exp 

P( {/?,.}) x { exp 
exp 



ft~2 Sr $r 

ft Sr $r ~ JjJ S<r,r'> S r S r '^j < fie < U . 



(14) 







Sr S r 



U < u c 



This distribution is a classical Ising model with a magnetic field term h J2 r Sr an d an 
antiferromagnetic spin-spin interaction — JJ2< r ,r'> S r S r > with J > 0. The latter is absent 
in the regimes with \x c < and /i c > U. These regimes are only controlled by an effective 
magnetic field /x//2 and (fif — U)/2, respectively: A positive field yields a positive Ising 
spin, implying that the optical lattice is fully occupied by heavy atoms. A negative field 
yields a negative spin, implying the absence of heavy atoms in the optical lattice. However, 
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there is no real competition between the two atomic species. In the intermediate regime 
< n c < U, on the other hand, there is competition due to the appearence of the linear 
and the bilinear spin terms: 

Hl=^^Y,Sr-^J E S r S r , (15) 
r <r,r'> 

The linear term favors a homogeneous distribution of heavy atoms whereas the antiferro- 
magnetic bilinear term favors a staggered distribution of heavy atoms (i.e. a site that is 
occupied by a heavy atom has nearest neighbor sites without heavy atoms). The latter is 
caused by the tunneling of the light atoms, since it is proportional to the square of tunnel- 
ing rate of the light atoms i. A simple mean-field calculation reveals that there is a first- 
order phase transition between the homogeneous phase and the staggered phase. Thermal 
fluctuations lead to the competition of homogeneous and staggered clusters A typ- 

ical realization of the distribution of heavy atoms, created by a Monte-Carlo simulation, 
is shown in Fig. 1. 

At \if = n c = U/2 the distribution in Eq. (JHJ) has a global Ising-spin flip symmetry and 
is related to a half-filled system. This means that half of the sites of the optical lattice 
are occupied by heavy atoms. In this case there is a continuous phase transition from 
antiferro- to paramagnetic order if the temperature is increased. 



5. Discussion 



The complex interplay between light and heavy fermionic atoms in mixtures results 
in a correlated distribution of the heavy atoms, given by the distribution density of Eq. 
(jHJ). The correlations are a consequence of the tunneling processes of the light atoms, 
i.e. they are due to quantum effects. If quantum tunneling is suppressed (e.g., by high 
barriers in the optical lattice and/or by using heavy atoms), the distribution of the heavy 
atoms looses the correlations, as discussed in Sect. 4.1. This is a consequence of our 
simple model where only local interatomic interactions are considered. In the presence 
of a nonlocal interaction between the atoms correlations would also exist in the absence 
of quantum tunneling. However, it is unlikely that nonlocal interactions play a crucial 
role for neutral atomic mixtures in an optical lattice because it is dominated by s-wave 
scattering. 

The correlations of the local density fluctuations of the heavy atoms lead to different 
phases and phase transitions. A perturbation expansion for weak tunneling (i.e. tunneling 
t is weak in comparison with the interatomic interaction U), presented in Sect. 4.2, has 
revealed that the heavy atoms are homogeneously distributed with one atom per optical 
lattice site if their chemical potentials fi c and fif are larger than the interatomic interaction 
U . According to the Ising Hamiltonian Hi of Eq. (|15|) . which describes the distribution 
for < fA c < U, a homogeneous distribution of heavy atoms also exists if fif — fi c is large 
in comparison with P/AU. On the other hand, for \fif — fi c \ small in comparison with 
P/4U the heavy atoms are arranged in a staggered order, similar to a charge-density 
wave in solid-state physics [T3j with rif ~ 0.5. For fij — ji c <C —P/AU heavy atoms are 
pushed out of the optical lattice. These results indicate that quantum tunneling causes an 
effective repulsive nonlocal interaction between the heavy atoms. The Ising Hamiltonian 
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Hj yields a first order transition from the homogeneous to the staggered distribution for 
a decreasing fij — ji c > 0. The associated coexistence of clusters with both types of order 
(see also Fig. 1) is similar to phase separation, discussed in the solid-state literature jl!2j . 
Furthermore, there is a second order transition from the staggered distribution at low 
temperatures to a disordered distribution at higher temperatures. In experiments the 
phase transitions are easily accessible by a variation of the tunneling rate i in the optical 
lattice. 

The symmetric situation \i$ = fi c , where the optical lattice does not distinguish between 
the heavy and the light atoms energetically, there is an additional type of competition 
between heavy and light atoms even in the absence of quantum tunneling: although there 
is no order in the distribution of heavy atoms, the sites of the optical lattice can be 
randomly occupied either by light or by heavy atoms, provided the repulsion is strong 
enough (i.e. for < \i c = fj,f < U). Consequently, the average density is n c = rif = 0.5 (cf. 
Eq. (fT3*|) ). The total symmetric case /if = fi c and if = i c is the Hubbard model. At half- 
filling the strong-coupling expansion leads to an antiferromagnetic spin-1/2 Heisenberg 
model [T3]. The ground state of this model is a staggered (Neel) state. If the spin is 
associated with two atomic species, the ground state of the atomic system is an alternating 
arrangement of the two type of atoms in the optical lattice. 

Light atoms are also affected by the order of the heavy atoms in the mixture, since 
they are scattered by the ensemble of heavy atoms. As a result, they can propagate (if 
heavy atoms are ordered, e.g., at low temperatures), they can diffuse (if heavy atoms are 
weakly disordered), or they are localized (if heavy atoms are strongly disordered). Due to 
correlations of the distribution of heavy atoms, a gap in the density of states of the light 
atoms can be opened |12|13| . This leads to an incompressible state of light atoms, where 
diffusion is absent. For disorder of heavy atoms (i.e. in the regime of competing clusters 
or at higher temperatures) localization of light atoms can take place. 

6. Conclusions 

There is a complex interplay between light and heavy fermionic atoms, where the latter 
can be described by classical Ising spin. The nonlocal interaction between the heavy 
atoms is caused by quantum tunneling of the light atoms. This is given by an effective 
Ising model with a symmetry-breaking magnetic field and an antiferromagnetic spin-spin 
coupling, at least in the regime of strong coupling between the two atomic species. It 
implies a competition between homogeneous and staggered order in the distribution of 
heavy atoms. There is a first order transition from homogeneous to staggered order and a 
second order transition from staggered order to a disordered distribution of heavy atoms. 

Acknowledgement: I am grateful to S. Fialko for providing Fig. 1. 
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Appendix A 



Using v = fi c — U/2 and g = U/2 we can write the determinant in equation 
power series of i (i.e. in the strong-coupling regime): 



cLS cL 



exp 



{Tr [in (l 



j(u+i- g s) 



} « det (Gr/ 1 ) exp 



Tr(dD) - -Tr(GiDGiD) 



with Gj l = 1 + e^ u ~ 9S) and D = e p{u+i - gS) - e 



A lenghty but straightforward calculation gives with A r = v — gS r the relation 



1„ f3^eP A * 



t 1 



TrfGxD) - -TrfGxDGiD) « - Y .. .. . 

y ' 2 y 1 1 2ff A r -A r , (1 + ^0(1 + ^) 



Since there are only values S r 

e f3A r _ e j3A r , 



E(^S r , S r 



-1, 1, the term 
P 



A r -A r > (1 + e^)( 1 + e' 3A r') 
can also be expressed as a quadratic form with respect to the Ising spins: 
E(S r , S r ') = Eq + E±(S r + S r r) + E2S r S r i 
with 



E = ~[£(l,l) + £(-l,-l)+2£(l,-l)] 



t 2 



'P 2 



+ 



(1 + 6^-/90)2 (i + efb+to) 



+ 







e pv sinh(/fy) 



Ag 1 + e 2 ^ + 2e^ cosh(/?^) 



£ 2 = -[£(l,l) + £(-l,-l)-2£(l,-l)] 



F 2 



Pv~Pg 



,Pv+Pg 



+ 



(l + e^"^)2 (i + e ^ff) 



/3 



eP v sinh(/3#) 



4# 1 + e 2 ^ + 2e^ cosh(/fy) 
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Figure 1. Competition of homogeneous and staggered clusters of heavy fermionic atoms 
(indicated as dark squares): a typical realization from the distribution in Eq. (jSJ) at 
Vf = Vc = U/2. 



and 

1 3 2 P 
E 1 = l [E(l,l)-E(-l,-l)] = Z± 

A further contribution comes from 



det (GT 1 ) = II I 1 + ef3{U ~ gSr 



(I _|_ e Pv-Pgy (I _|_ e flv+0gy 



